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Experimental evidence of electron spin precession during travel through the phosphorus-doped Si 
channel of an all-electrical device simultaneously indicates two distinct processes: (i) short timescales 
(~50ps) due to purely conduction-band transport from injector to detector and (ii) long timescales 
(~lns) originating from delays associated with capture/re-emission in shallow impurity traps. The 
origin of this phenomenon, examined via temperature, voltage, and electron density dependence 
measurements, is established by means of comparison to a numerical model and is shown to reveal 
the participation of metastable excited states in the phosphorus impurity spectrum. This work 
therefore demonstrates the potential to make the study of macroscopic spin transport relevant 
to the quantum regime of individual spin interactions with impurities as envisioned for quantum 
information applications. 



Incorporating impurities into the otherwise pure silicon 
(Si) lattice has numerous consequences for charge trans- 
port. Most importantly, doping using atomic species with 
more than Si's four valence electrons not only provides 
mobile charge to the conduction band, but also leaves 
behind a positively-charged ion that modifies the elec- 
trostatic energy landscape when not screened. But be- 
yond impacting the flow of electron charge, electron (or 
n-type) doping also impacts spin-polarized transport as 
well. 09 

For example, several years ago we showed how band 
bending gives rise to non-ohmic spin transport in «3 
/im-thick n-type lightly phosphorus-doped Si using all- 
electrical ballistic hot electron injection and detection 
techniques. [H In the present work, we show that un- 
der certain circumstances, long-timescale processes not 
seen in otherwise-equivalent undoped devices [jj @ can 
be observed in measurements of devices with lowly-doped 
transport channel regions. These results are attributed 
to the interaction of conduction electrons with shallow 
impurity-related traps and suggest that the study of spin- 
polarized transport in semiconductors can potentially be 
used to elucidate physics previously accessible only to 
time-domain techniques or to explore the local interac- 
tion of spin information with isolated impurity potentials. 

The devices we use here are nominally identical to 
those in Ref. 0, but with slightly lower impurity density. 
As illustrated in Fig. [1] (a) , spin-polarized hot electrons 
are injected from a ferromagnetic (FM) tunnel-junction 
cathode (CoFe, biased at emitter voltage Ve) and travel 
ballistically through a nonmagnetic thin film and over 
an energy-filtering Schottky barrier to couple with con- 
duction band states in the 3-micron-thick n-Si transport 
layer. Spin detection after vertical transport through this 
layer (biased at collector voltage Vex) is accomplished 
by analyzing the ballistic component of the hot-electron 
current generated by ejection from the n-Si conduction 
band over a Schottky barrier and into a second FM thin 
film (NiFe) . Because the mean- free-path of hot electrons 
is determined by the relative orientation between their 
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FIG. 1: (a) Schematic band diagram of vertical n-type doped 
Si spin transport channel devices used in this work, show- 
ing band-bending and the resultant potential minimum for 
temperatures high enough to ionize dopants and low applied 
voltage Vci- (b) Illustration of model implied by Eqn. [T]for 
capture of electrons (incident with period n) by N shallow 
traps AE below the conduction band with probability a. Each 
trap is initially empty with probability Pe and releases the 
electron with an exponential probability distribution having 
timescale Tt. 



spin and the magnetization of the FM film, the magni- 
tude of this current (Ic2) then carries information about 
spin transport in the channel. 

Previous measurements on these devices have indi- 
cated that the confluence of Schottky depletion regions 
on both injector and detector sides of the transport chan- 
nel result in a confining conduction-band profile (i.e. an 
electrostatic potential energy minimum exists) for ap- 
plied voltages between injector and detector (Vci) be- 
low s»2 V at temperatures > 40 K, sufficient to ionize 
the dopants. Assuming full depletion, this biasing be- 
havior corresponds to a doping density of approximately 
3xl0 14 cm^ 3 . 
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FIG. 2: Perpendicular-field measurements at 55K. In (a), a 
full field loop displays square hysteresis and clear low-field 
dephasing. In (b), voltage dependence of the low- field signal 
is shown, where the applied voltage V c i is changed from OV, 
0.5V, 1.0V, and 1.5V to 2.5V in steps of 0.1V. 

In a magnetic field oriented perpendicular to the de- 
vice plane, the magnetization of injector and detec- 
tor magnetic thin films remains in-plane due to shape 
anisotropy. However, spins injected into the nonmag- 
netic Si are induced to precess around this magnetic 
field; in drift-dominated transport through the conduc- 
tion band, the transit time is well-defined and this leads 
to coherent spin-signal oscillations as a function of mag- 
netic field magnitude, jf| H, Q On the other hand, when 
a large transit-time uncertainty exists (due e.g. to 
random diffusion), electron spins arrive at the detec- 
tor with a distribution of precession angles. This spin 
"dephasing" suppresses the oscillation amplitude. A 
Fourier-transform method applied to the quasi-statically 
measured magnetic-field spectroscopy can directly re- 
cover the empirical spin transit-time distribution with ns 
resolution (jj fl(]| and indicates an inverse relationship be- 
tween the magnitude of magnetic-field features and trans- 
port timescales. 

Measurements on our phosphorus-doped Si devices 
with a perpendicular magnetic field (at T=55K in Fig. [5] 
(a)) show an open loop corresponding to in-plane mag- 
netization switching of the softer NiFe detector layer (in- 
duced by small field misalignments [llj]), superimposed 
on the signal due to spin precession during transport. 
It is evident that our spin transport signal includes the 
signatures of two distinct processes: both (i) precession 
features on the kOe scale beyond the range shown due to 
the short («50ps) timescale associated with conduction- 
band transport from injector to detector, and (ii) a low- 



FIG. 3: Temperature dependence of low-field dephasing be- 
tween 40K and 100K at applied voltage Vci =1.5V. (b) 
Simulated spin-precession Hanle signals using the trapping- 
time model (Eq. [1} with parameters a— 0.2, ro=200ps, 
A_E=11.5meV, and effective 1-d current of 3nA sampling 
N =50 traps. Inset: Calculated trapping-time distributions. 

field dephasing resulting in a narrow peak of width 50- 
100 Oe due to much longer timescales. This observation 
is further supported by noticing that the small-field peak 
in a parallel injector-detector magnetization configura- 
tion (left-to-right sweep, red) becomes inverted for the 
antiparallel configuration (right-to-left sweep, blue). 

This small-field peak is sensitive to the applied voltage 
between injector and detector (Vci)- In Fig- IHb), we 
show the low-field region of our measurements from par- 
allel magnetization (Jp), where changes in absolute signal 
current as a function of applied voltage are eliminated by 
plotting spin polarization (Ip — I a vg) / lavg, where I aV g is 
the average between parallel and anti-parallel signals, as 
shown by the horizontal dotted line in Fig. HJa). Clearly, 
higher applied voltage corresponds to a widening peak 
that accounts for less of the total spin polarization sig- 
nal. At voltages sufficient to entirely remove the confin- 
ing potential in the conduction band (« Vci =2.5V), the 
low-field peak disappears entirely. 

Temperature dependence of this phenomenon at a fixed 
Vci =1.5V is shown in Fig. Efa). Here, we see that 
increasing temperature reduces the relative contribution 
of the low-field peak to the overall signal and obscures 
the simultaneous presence of large-field scales. Obser- 
vation of the dual-timescale phenomena it otherwise im- 
plies therefore requires a narrow range of both voltages 
and temperatures; this explains why it was overlooked 
during the study which led to Ref. |j. 

This low-field feature is reminiscent of similar obser- 
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vations in optical[l2| and transport[l3} measurements 
of spin-polarized electron precession in GaAs. In that 
case, the phenomena was attributed to a non-linear dy- 
namic feedback between electron spin and nuclear spin 



via Overhauser and Knight fields. [14j Ref. [l3j, for in- 
stance, demonstrates a remarkable degree of correspon- 
dence with steady-state finite-difference time-domain cal- 
culations of spin transport in the presence of the self- 
consistent nuclear field, and its modification by mi- 
crowave excitation resonant with the nuclear spin split- 
ting. 

Despite the similarity to the results observed here 
with Si, however, there are several reasons to disregard 
this hyperfine-mediated mechanism as a possible origin. 
Si is expected to have much weaker hyperfine interac- 
tions than GaAs due to the small relative abundance 
of nonzero nuclear spin isotopes (~5% for spin-1/2 29 Si; 
spin-1/2 31 P dopants are even more dilute). In compar- 
ison, every nucleus in GaAs has nonzero spin. Further- 
more, we do not observe the measurement time depen- 
dence associated with long nuclear spin lifetimes seen in 
GaAs spin transport devices. [iH ITBIj 

Our spin current signal is essentially a measurement 
of the average spin precession angle 0, which is a prod- 
uct of both spin precession frequency uj and transit time 
t. In principle, the effects discussed above could be the 
result of modification of either of these; a linear process 
would require that there is a sub-ensemble of electrons 
either moving substantially slower or precessing substan- 
tially faster than the rest to account for both large- and 
small- field features. The latter cause is eliminated by the 
nonexistence of g-factors significantly greater than 2 in 
Si. 

The origin of our observations is therefore a transit- 
time effect. Here, we establish that this is due to the 
presence of trapping into and subsequent re-emission of 
spin-polarized conduction electrons from shallow phos- 
phorus impurity-related states. At typical injection cur- 
rents of «10 fj,A, there are approximately 10 14 electrons 
entering (and leaving) the transport channel each sec- 
ond. For a quasi-Lorentzian Hanlc half-width at half- 
maximum of AB «50 Oe, an average transit time of 
h/ (<7/isAS) wl ns implies a steady-state population of 
10 5 occupied traps in the channel. This value is below 
the actual density of approximately 10 7 available charged 
impurities under conditions above the ionization temper- 
ature in the effective volume «100/^mx 100/mix3/mi at 
the phosphorus density used here (although this num- 
ber must be reduced at low temperatures due to non- 
negligible occupation probability). 

By varying the electron density via injected current 
(Ici)i we can see the signature of this fixed trap num- 
ber. Although the available range of this parameter is 
limited by the reliable operating conditions of our tunnel 
junction spin injector, one can clearly see in Fig. HJa) 
that the overall polarization increases with injected cur- 
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FIG. 4: Low-field dephasing as a function of injected current 
and hence steady-state electron density. In (a), spin polar- 
ization is shown, whereas (b) demonstrates a reduction in 
relative contribution to the overall signal from trapped elec- 
trons when density increases. Our numerical model included 
for comparison uses an effective 1-dimensional current. 



rent. The relative contribution of the low- field peak from 
this data is shown in Fig. E^b), from which it is appar- 
ent that a small but substantial decrease in the relative 
proportion of trapped electrons (i.e. relative peak height 
to background polarization) accompanies an increase in 
conduction electrons as well. This behavior indicates the 
approach to saturation of the impurities and a fixed con- 
tribution to low-field dephasing. 

Average transit times for electrons participating in 
trapping events of approximately Ins (indicated by the 
width of the low-field peak) are relatively small at these 
temperatures for the expected ground-state phosphorus- 
impurity trap depth of 45meV.[l6| Consistency of the 
present trapping scheme with the observed timescales re- 
quires the participation of trapping-emission cycles in- 
volving shallow states that are presumably excited Ry- 
dberg states with suppressed relaxation to the rela- 
tively deep ground state. In this regard, it is sup- 
portive to note that recent time-domain experiments 
with a far-infrared pulsed free-electron laser has iden- 
tified transition lifetimes of the phosphorus 2p excited 
state (AE =11.5meV below the conduction band) ex- 
ceeding 200ps in this temperature range. [I?} Further- 
more, recent calculations su gges t an intrinsic 2p-ls tran- 
sition timescale of over Ins, 18] which has even enabled 
population inversion and lasing in externally-pumped 
systems. 

We now describe a simple numerical model to calculate 
the transit-time probability function which incorporates 
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the details of trapping into (and emission out of) the 
impurity state. To capture the observed behavior, one 
must consider the convolution of all < k < N possi- 
ble trapping events in N available traps with trapping 
probability a: 



N 

■E 

fe=i 



(1 - aP E ) N S(t) + V ) (1 - aP E ) N - k (aP E ) k G{t, k, r t ) 



(1) 

where G(t,k,Tt) is the gamma probability distribution 
(convolution of k exponentials of trapping timescale Tt). 
The probability of an empty trap (Pe) is given by the rel- 
ative ratio of incident timescale (ji determined by the in- 



jection current), to trapping timescale (r f = roe 



AE/k B T 



such that P E = 1 



The physical significance of 



Ti/a+r t ' 

the relevant parameters of this model are illustrated in 
Fig. ffl(b). 

Since this a one-dimensional model, an equivalent "ef- 
fective" current (to determine n) and trap number (iV) 
must be used. At the phosphorus density used here, the 
inter-dopant distance is approximately 150nm. There- 
fore, with a total electron injection area of 10 4 /zm 2 , 
w 2.25 x 10~ 6 of the injected current will interact with 
N « 20 dopants through the 3-/im transit length. A real 
injected current Ici — 79/iA at Ve = —1.5V is therefore 
equivalent to an effective 1-d current of approximately 
200pA. Allowing for some degree of variation in these 
parameters given the simplicity of the analogy, we use an 
effective current of 3nA and N = 50 in our subsequent 
comparison of the model fit to the experimental data. 

Ignoring the contribution from the much smaller 
timescale features endowed by drift/diffusion/relaxation, 
the observed spin signal is then the real part of the 
Fourier transform of the distribution in Eqn. [TJ0, Hoj 
Correspondence of this model's numerical predictions to 
the general features of experimental data are evident in 
Fig. [3]Jb), where the changes in trapping contribution 
and weak linewidth dependence on temperature observed 
in experimental results in Fig. [3][a) are clearly repro- 
duced. Furthermore, we also show that the numerical 
model captures the dependence of the relative contribu- 
tion from trapped electrons on injected current as seen in 
Fig. H] (b) using identical model parameters. Although 
we have not included it in the model, the broadening and 
peak reduction seen in voltage dependence (Fig. [2]) can 
be attributed to field emission and a subsequent reduc- 
tion in trapping timescale r t and capture probability a. 

The explicit coupling of spin-polarized conduction elec- 
trons with phosphorus impurities as demonstrated here 
opens many research possibilities. It may be useful 
as a probe of impurity levels and transition rates not 
explicitly requiring time-domain methods as previously 
assumed were necessary. 17( In addition, we speculate 
that it may become possible to explore the role of 
contact hyperfine interactions and dynamic nuclear po- 



larization of the phosphorus nuclear spin by itinerant 
non-equilibrium spin-polarized conduction electrons with 
NMR techniques forpotential application to quantum 
computing schemes. [201] 
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